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^ Abstract. The classical continued fraction is generalized for studying 

the rational approximation problem on multi-formal Laurent series in 
this paper, the construction is called m-continued fraction. It is proved 
that the approximants of an m-continued fraction converge to a multi- 
formal Laurent series, and are best rational approximations to it; con- 
versely for any multi-formal Laurent series an algorithm called m-CF 
transform is introduced to obtain its m-continued fraction expansions; 
' moreover, strict m-continued fractions, which are m-continued fractions 

imposed with some additional conditions, and multi-formal Laurent se- 
ries are in 1-1 correspondence. It is shown that m-continued fractions 
can be used to study the multi-sequence synthesis problem. 



^ . 1. Introduction 

I— I . 

' Continued fraction ^J[^|211llHl[ini is a useful tool in dealing with many 

number theoretic problems and in numerical computing. It is well-known 
I that the simple continued fraction expansion of a single real number gives 

■ the best solution to its rational approximation problem. Many people have 

contrived to construct multidimensional continued fractions in dealing with 
i -^H I the rational approximation problem for multi-reals. The archetypal exam- 

d • pie of a multidimensional generalization of simple continued fraction is the 

Jacobi- Perron algorithm (JPA), see [2]. This algorithm and its modification 
are extensively studied ^1 El El El • These algorithms are borrowed to 
study the same problem for multi- formal Laurent series [HI • But none 
of these algorithms guarantee best rational approximation in general. In 



X 



■ this paper, we deal with the multi-rational approximation problem over the 

formal Laurent series field F{{z~^)): given an element r G F{{z~^))^ , find 
p G € F[z], such that p/q approximates r as close as possible while 

deg{q) is bounded. In this setting, we give a natural generalization of the 
simple 1-dimensional continued fraction to the multidimensional case, and 
demonstrate that approximants of the continued fraction expansions of a 
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multi-formal Laurent series r we define indeed give the best rational esti- 
mates of r. 

The rational approximation problem of the multi-formal Laurent series is 
motivated by sequence synthesis problem, which has applications in the field 
of communication and cryptography. It is known that the single-sequence 
synthesis problem can be solved by the iterative Berlekamp-Massey (BM) 
algorithm ^ j the Mills algorithm using continued fractions [T^] , and the 
Euclidean algorithm as presented by Sugiyama Y. et al. [22j. Some con- 
sequential work shows that the continued fraction algorithm is a powerful 
tool [SI El 1201 HI] • Though for solving sequence synthesis problem the con- 
tinued fraction algorithm is equivalent to the BM algorithm and Euclidean 
algorithm, some mysterious data structures of the latter have natural inter- 
pretation in view of the former [IJ [3] . The BM algorithm and the Euclidean 
algorithm can be generalized to solve the multi-sequence synthesis prob- 
lem i6j|7j, which can also be solved by the lattice basis reduction algorithm 
as presented by Wang L.P. et al. HHl HHl • It is natural to think of the 
generalization of continued fractions to multidimensional case. 

Our main contributions in this paper are: 1, we give a definition of contin- 
ued fraction in the multidimensional case; 2, set up the transformation from 
continued fractions to elements of -F((z~^))'"; 3, prove that the approximants 
of a continued fraction are best rational approximants of its corresponding 
element; 4, give the procedures of computing the continued fraction expan- 
sions of a given element of F((z~^))™; 5, give a natural 1-1 correspondence 
between strict m-continued fractions, which are m-continued fractions im- 
posed with some additional conditions, and elements of F{{z~^))'^ , 6. apply 
the theory to solve the multi-sequence synthesis problem. 

The rest of this paper is arranged as follows: Section 2 lists the notations 
used in the paper. Section 3 deals with the indexed valuation of F{{z~^))'^. 
Section 4 is the detailed definition of the problem of multidimensional ra- 
tional approximation. Section 5 defines m-pre-continued fractions and m- 
continued fractions and gives some of their basic properties. Section 6 states 
the main results of this paper. Section 7 completes all proofs. 



F denotes a field, F{{z ^)) = <^Ylt>i '^t^ * | ^ G .Z', at G F | is the formal 
Laurent series field [H], where Z denotes the integer ring. F[z) denotes 
the rational fraction field over F, F[z\ denotes the polynomial ring over F. 
The discrete valuation on F{{z^^)) is denoted by v{-), i.e., v{J2t>i ^t^^"*) = 
i if a, / 0. For a = Ed>-t>o a*^"* + Ei<t<oo«t^"* ^ F2{{z~^)), the 
polynomial part X]^>_t>o is denoted by [aj ; and the remaining part 

Z]i<t<oo "-t^^^ is denoted by {a}. 

m denotes a positive integer, Zm denotes the finite set {1,2,- •• ,m}. 
F{{z~^))^ denotes the set of all m-tuples over F{{z~^)), similar for F{z)"^ 



2. Notations 




MULTIDIMENSIONAL CONTINUED FRACTION AND RATIONAL APPROXIMATIONS 



and ^[z]™, and we always write those m-tuples in column form. For an ele- 
ment r = {ri{z),r2{z), • • • , rm,{z)y in F{{z~^))"^ , where r means transpose, 
we denote {lri{z)\, lr2iz)\,--- , km(^)J)^ by [rj, and 
{{niz)},{r2{z)},--- ArmizW hj {r}. 

In denotes the identity matrix of order n for any positive integer n, and 
Cj denotes the j-th column of the matrix /„ for j : 1 < j < n. Denote by 
Ma^biFiiz~^))) the set of all possible matrices of order a x b over F{{z~^)) 
for any integers a and b, and similar for M(j fe(F[z]). 

3. Indexed valuation on F{{z~^))^ 

In this section the indexed valuation on F{{z^^))"^ is introduced, which 
will be in use in studying the rational approximation problem for multi- 
Laurent series. 

Definition 3.1. ( Linear ordering on x Z ) 

For any two elements {h,v) and {h',v') in Zm x Z, we say {h,v) < {h',v') if 
V < v' , oi V = v' and h < h' . 

It is clear that the ordering defined above on Z^a x Z is linear ^lOj,, i.e., 
any two elements can be compared. 

Lemma 3.2. Let {i,x) and {j,y) be any two elements in Z^i x Z, then 
{i, x) < (j, y) if and only if x + Lij < y, or equivalently if and only if 
x-ljA < y, where kj = 1 ifi>j, hj = Oifi< j; Lij = lifi>j, Lij = 
ifi < j- 

Proof. It is easy to check. □ 

For any {j,x) G Z^ x Z, we call the non-zero element z^ej in F{{z~^))'^ 
a monomial, and denote 

(3.1) Iviz''ej) = {j,-x) eZmX Z. 

Notation 3.3. Denote by M the set of all monomials in F{{z~^))"^. Any 
element r in F{{z~^))"^ can be expressed uniquely as a sum of the following 
form: 

(3.2) r = ^ Cmni,Crn^F,m^M. 

We will say m is a monomial of r and denote m € r if Cm 7^ 0. The total 
number ( finite or infinite ) of the nonzero coefficients Cm will be called the 
Hamming weight of r, denoted as W}{{r) = t- 

Definition 3.4. {Indexed valuation) 

For an nonzero element r = {ri(z),r2{z), ■ ■ ■ ,rm{z)Y in F{{z~^))'^ , we call 
Iv{r) = min{/T;(m) | m G A4, ?tt, G r } G Zm x Z 
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the indexed valuation of r, where Iv{rn) is defined as in If Iv{r) = 

{h,v), then v will be called the valuation of r and denoted by v{r), and /i 
will be called the index of r and denoted by /(r). It is clear that 

/ ^(n) = min{f (rj(z)), 1 < i < m}, 

\ ^{l) — min{i[l <i< m,v{ri{z)) = v{r)}, 

where v{ri{z)) is the discrete valuation on F{{z~^)) 9 . By convention 
v{0) = oo, lv{0) = (1, oo), and Iv{r) < lv{0) Mr / 0. 

For r € F{{z^^))'^ such that wh{z) = t < oo, define 

Supp{r) = {Iv{m)\rn G A4,rn, € r}, 

which will be called the support of r; and denote by Supp~^{r) the largest 
element in Supp{r) if wh{l) < The following lemma is clear. 

Lemma 3.5. The indexed valuation is a surjective mapping Iv: F{{z~^))"^ 
{Zm X Z) L) (1, oo) which satisfies 

(1) Iv{r) = (l,oo) <^ 11= 0. 

(2) Iv{ru{z)) = {h,v + v{u{z))) for all r G Fiiz-^))"" and u{z) G 
F{(z~^)), where {h,v) = Iv{r). 

(3) Iv{r + s) > mm{Iv{r), Iv{s)} for allr,s G F{{z^^))^ . The equality 
holds true if Iv{r) ^ Iv{s). 



4. Rational Approximation 

In this section we give a detailed definition for the rational approximation 
problem on multi-formal Laurent series. Let r G F{{z~^))"^, ^ q{z) G 

F[z], p{z) G Flz]'"^, we call a rational approximant to r if v{r — "^(^) > 

deg(g(z)), and call Iv{r — ^(^) the precision of to r, here it is worth 

of pointing out that the components of have a common denominator 

q{z). It is easy to see that is a rational approximant to r if and only 
if p{z) = [rq{z)\ . So any nonzero polynomial q{z) can be the denominator 
of a rational approximant to r, and such a rational approximant is unique. 
As a matter of convenience, any nonzero polynomial q{z) will be called a 
denominator of r with precision Iv{r — ^^||fj~^)- 

Definition 4.1. {Best Rational approximant) 

Let r G Fffz"^))'" . We call $1 € F(z)"^ a 6esi rational approximant of r 
if denominators of r with degree lower than deg{q{z)) have precision lower 
than Iv{r_ — ^py), and denominators of r with degree same as q{z) have 



denominator of r. 



precision no greater than Iv{r — =||y)- In this case, q{z) is also called a 6esi 
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Let S be the set of all elements r of F{{z ^))'^ with v{r) > 0. For 
the rational approximation problem on F{{z~^))"^ , we need only consider 

elements in S. In fact, if we write r = [rJ then, {r} E S, and [rJ +|||y 

is a rational approximant to r of precision (h, n), if and only if =||y G F{z)"^ 

is a rational approximant to {r} of the same precision {h,n); and q{z) is a 
best denominator of r if and only if it is a best denominator of {r}. 

In studying the best rational approximation problem for r € S, motivated 
by the simple continued fractions in the case m = 1, we will introduce the 
multidimensional continued fraction in the next section. 

5. Multidimensional continued fractions 

In this section we introduce the m-pre-continued fractions and the m- 
continued fractions, and give some of their basic properties. 

5.1. m-pre-continued fractions. 
Definition 5.1. {m-pre-continued fraction) Let 

(5.1) C = [ao,/ii,ai,/i2,a2, • • • , /ifc,afci • ' 'l^ < A; < cj, 

where G and Of. G ^[z]™ for all /c : < A; < u, and w is a positive 
integer or oo, then C will be called an multidimensional-pre-continued frac- 
tion, or simply m-pre-continued fraction, and iv the length of C. We always 
assume ag(z) = 0, since we will see ag is irrelevant to our concern. 

Notation 5.2. Denote by the matrix of order (m + 1) which comes 
by exchanging the h-th column and the (m + l)-th column of the identity 
matrix Im+i- 

(5.2) Eh = (ei 62 • --e^.i e.^+i e^+i ■ ■ ■ C/, ) G Mm+i,m+iiF[z]). 

From C defined as (|5.1|1 we define iteratively the square matrices of 
order (m + 1) over F[z]: 

(5.3) Bo = Im+i, Bk = Bk-iEh^A{ak) G Mm+i,m+i{F[z\), k>l, 
where 

(5.4) A{ak) = ( -^'^ ) G M„+i,^+i(F[z]); 

and denote by by. the last column of B^, and let p^ = Pj^{z) G F[z]"^ and 
Qk = Qkiz) G F[z] be components of 6^, that is 

(5-5) h = iPk^QkY- 

Some properties of the sequence {{Pj^, qk)}k>o are given in the following 
proposition. 

Proposition 5.3. (1) gcd{qk{z),pk,i{z),pk,2{z), ■ ■ ■ ,Pk,m{z)) = 1 for /c > 
0, where Pkjiz) is the j-th component of p^. 
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(2) For A; > 0, let Pfc-i be the matrix of size mxm and Qk-i the matrix 
of size 1 X m such that 

(5.6) B>^=(o'' ? 

and denote by Pj!j_]^j(€ F[z]™) the jth column of and Qk-ij{& 
F[z]) the jth component of Qk-i, 1 < i < "i. Then for k > 1, 

qk J V Qk-i J-""- ' \ Qk-2,hk 
or explicitly 

Ik = Pk^i"'k,hkiz) + Yl Pk-2,jak,jiz) + Pk_2,hk, 

jV/ifc,l<j<m 

Qk = qk-iak,hki^) + X] Qk-2,jak,j{z) +Qk-2,hk- 

jjthk,l<j<m 

(3) For 1 < k < to and j S Zm, denote by l{k,j) the largest positive 
integer k' < k such that h^i = j, and let l{k,j) = if no such k' 
exists. Then for k > 1, 

Proof. It is easy to prove. □ 
5.2. Conditions 1-4 and m-continued fractions. 

Definition 5.4. {Conditions 1-4) For the m-pre-continued fraction C de- 
fined as (|5.1|) . we define condition 1 as below: 

• Condition 1: deg{ak,hk{^)) > 1, V 1 < < w, where ak^h^{z) denotes 
the /ifc-component of aki^)- 
In the sequel, we always assume C satisfies the condition 1, and associate it 
with the following quantities (for each A; : 1 < A; < w): 

{ifc = deg(afc,/j^(z)), to = 0, 
'^kj = Eh,=j,i<k'ti^ '^o,j=0, 
Vk=Vk,hk, vq = Q, 

and the diagonal matrix Dk (for each A; : 1 < A; < cj): 

(5.8) Dk = Diag.{z- 



\ 



We define conditions 2-4 on C as below (by convention oo — 1 = oo): 

• Condition 2: (/ifc, ffc-i.^J < (/ifc+i, f^+i), 1 < A; < u; - 1. 

• Condition 3: Iv{Dkak) = ihk,Vk-i^hk)i 1 < A; < w. 
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• Condition 4'- Supp^{Dkaf.) < {hk+i,Vk+i), 1 < k < uj — 1. 

Definition 5.5. An m-pre-continued fraction C is called an m-continued 
fraction if C satisfies the conditions 1-3. An m-continued fraction C is said 
to be strict if it satisfies the condition 4. 

It is clear that Conditions 1, 3 and 4 imply Condition 2. These conditions 
can be stated in some equivalent forms as shown in the following proposition. 

Proposition 5.6. Let C, which is defined as (|5.1j) . satisfies the condition 
1, keep all notations made for it, and denote 

(5.9) = {ak,i{z), ■ ■ ■ ,afcj(z), • ■ ■ ,ak,m{z)y, I <k <uj, 

where akj{z) G F[z]. Moreover, we associate C with some quantities as 
below: for each k : 1 < k < co, 

( do = 0, 

(5.10) <^ dk = Ei<i<ktu 

[ nk = 4-1 + Vk{= 4 + Vk-i,h^:), 

where the notations and Vk are defined in 1)5. 7|) . Let and Ljj be defined 
as Lemma 13.21 Then 

(1) For k : 1 < k < UJ — 1, the following conditions are equivalent: 

(a) {hk,Vk-i,hk) < {hk+i,Vk+i)- 

(b) {hk,nk) < (/ifc+i,nfc+i). 

(c) Vk-l,hk ~ '^k,hk+i + ^hk,hk+l ^ 

(2) For k : 1 < k < UJ, the following two conditions are equivalent: 

(a) IviDkOf,) = {hk,Vk-i,hk). 

(b) deg(afc,j(z)) < Vkj - Vk-i,hk - Ihkj, VI < j < m,j / hk- 

(3) Assume the condition 3 holds true. For k : 1 < k < uj — 1, the 
following two conditions are equivalent: 

(a) Supp+{Dkaf,) < (/ifc+i, ffc+i). 

(b) For every j ■ 1 < j < m, akj{z) is of the form 

0<x<Xk,j 

where Xkj = mm{tk,j,Xkj}, tkj = Vk,j - Vk-i^h^ - h^j-, ^kj = 
Vk+i - Vk-i^hk - IhkJ - Lj,hk+i- 

(4) Assume both of the condition 2 and 3 hold true. For k : 1 < k < uj 
and j ^ hk, 

( deg{akj{z)) < dk- d/(A:,j)-i if K^J) > 1> 
\ deg{ak,j{z)) <0 ifl{k,j)=0. 

Proof. It is easy to prove. □ 

6. Main Results 
In this section we state the main results of this paper. 
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6.1. Approximants of m-continued fraction-I. In this subsection, we 
always let C, which is defined as (|5.1j) . be an m-continued fraction. We keep 
all notations made for it in section 5. 

Theorem 6.1. deg{qk{z)) = dk for all k : < k < u. 

Based on the above theorem, each pair (p^, qk) provides a rational fraction 

=^|^, which will be called the kth approximant (or, convergent) of the m- 
continued fraction C. The properties of these approximants are summarized 
in the following theorems. 




. ^ . X , ihk,nk), 1 <k <uj. 

qk-i{z) qk{z) J 

As a consequence, the sequence convergent in the case uj = oo. 



Denote 



if LU < OO, 



limfc^oo t^Vty, lo = oo. 



We see ^{C) € F{{z~^))'^, and call C an m-continued fraction expansion of 
f{C). Denote by C(r) the set of all possible m-continued fraction expansions 
of r . 

Corollary 6.3. Let r = f{C). For k : 1 < k < to, we have 

(1) Iv{r- = {hk,nk). As a consequence, rq^^iiz) - P,^_-^^{z) = 
{rqk-i{z)} and Iv{{rqk^i{z)}) = {hk,Vk). 

(2) Iv{r) = {hi,ti). 

(3) Iv{p^{z)) = (hi, —dk + ti). As a consequence, deg{p^{z)) = dk — ti, 
where deg{p^{z)) denotes the largest one among deg{pkj{z)), 1 < 

j < ITl- 

Proof. It is easy to prove. □ 

Theorem 6.4. Let r = (p{C). Assume q{z) G F[z\, dk < deg{q{z)) < dk+i 
for some < k < uj {d^j = oo if oj < oo). Then 

q[z) qk{z) 

As a consequence, each , < k < lo, is a best rational approximant 
to r, and the degree of any best denominator of r must be dk for some 
k : < k < UJ. 
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6.2. m-CF Transform and m-continued fractions-I. Given an element 
r S F[{z~^))"^ with r ^ and v{r) > 0. In this subsection we introduce 
a transform, which we call m-CF Transform, which may produce the m- 
continued fraction expansions of r. 

Definition 6.5. (D-matrix) We call a diagonal matrix over F{(z^^)) a D- 

matrix if each of its diagonal elements is a power of z. 

m-CF Transform: Given r G F{{z^'^))'^, £ 7^ and v{r) > 0. Ini- 
tially, take ao = G Flz'^]"", A 

— 1 — Im (the identity matrix of or- 
der m), and /3o = £• Suppose [aQ,hi,ai,h2,a2, ■ ■ ■ ; ^fc-i, afe-i], ^k~2 = 
Diag.{- ■ ■ , , • • • ) which is a Z?-matrix of order m, and = 

(■ ■ ■ • • • )^ ^ F((2:~-'^))™ have been defined for an integer k > 1. If 

/5fc-i = 0) the algorithm terminates and denote uj = k. If Pk-i 7^ 0) then do 
the following steps: 

(1) Take hk = /(Afc_2/3fc-i) S Z^. 

(2) Take Afc„i = Diag.{- • • , z'"'''^,- ■ ■ ), where c^j = Cfc_ij if j 7^ /i^, 
and Cfc,ft^ = v{Ak-2Pk-i) G ^• 

(3) Take = [/OfcJ -e^ and (3k = {pfcl+e^, where = (• • • , p^j, ■■■Ye 

F((z-i)r, p,,,- G p,,, = if J / p,,,, = 

Pk-ij is the j-th component of f3k-i, and G ^[z]™' is chosen freely, 

except that it satisfies the condition Iv{Ak-i{pk}) < -^^(^fc-iffc)- 

Denote tj = 00 if the above procedure never terminates. 

An m-CF transform on r will result in an m-pre-continued fraction: C = 

[0, /ii,a^, /i2, ■ ■ ■ ; ^fc; «/c5 ■ ■ ■ ]• Note that e^, may not be unique at each 

step k > 1, different choices of e^, will give different C. 

Definition 6.6. We denote by T(r) the set of all possible m-pre-continued 
fractions obtained from r by m-CF transforms. 

Definition 6.7. {A-polynomial part of r) Let r G F{{z~^))'^, A be a D- 
matrix. Write A = \r_\ = Ylm&M '^m IRj <^m & F, o = {r}, where A4 is 
defined in Notation 13.31 We call 

the A-polynomial part of r , where 

(6.1) 5nppA,a,-(^) = {meM\meA, Iv{Am) < Iv{Aa) }. 
Denote [zijj = CmlR, where 

meSuppA_a,+{A) 

(6.2) SuppA,a,+iA) = {rneM\meA, Iv{Aa) < Iv{Am) }. 

We see [rj = [rj^ + [nj^, by noting that the mapping m 1— > Am is a 
bijection on M.. 
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Lemma 6.8. Let r G F{{z-'^))"', A = [r\, a = {r}; and let A be a D- 
matrix. Assume e E F[z]"^, and a = A — e. Then the following conditions 
are equivalent: 

(1) «= LiJi- 

(2) Iv{Aa) < Iv{Ae) ifa = 0, and Supp+{Aa) < Iv{Aa) < Iv{Ae) if 
a / 0. 

{Supp{a) U Supp{e) = Supp{A), 
Supp{a) C Supp^^a-{A), 
Supp{e) C SuppA,a,+{A), 
where 5nppA,a,-(^) o-nd SuppA,a,+iA) are defined as \6. 1\) and i6.S\) . In 
particular, [rj^ = r for any r £ F[z]"^. 

Proof. It is easy to prove. □ 

Definition 6.9. {Mapping ip) Define ip{r) to be the element in T(r), ob- 
tained from r by the m-CF transform, where at each step we choose = 

L^fe^iSfc-2JAfe_i> we choose = L^fc^iSfc-2J A;,_i- 

The properties of m-CF transform are summarized in the following theo- 
rem. 

Theorem 6.10. Denote by S the set of all possible r € F{(z^^))'^ with 
v{r) > 0, and by C* the set of all possible strict m-continued fractions. 
Then 

(1) T(r) = C(r) VrG S. 

(2) The mapping ip is a bijection from S onto C* , and ip is its inverse. 

6.3. An application to multi-sequence synthesis problem. In this 
subsection we show m-continued fractions can be used to solve the multi- 
sequence synthesis problem. The multi-sequence synthesis problem was 
stated in terms of linear feedback shift registers [HI [7j. For convenience, 
we restate it by means of the indexed valuation on For any 

given sequence r = {ct}t>o over F, where q G F, we identify it with the 
formal Laurent series r(z) = X]t>o with valuation larger than 0, 
and let r^") = (co,ci,-- - ,c„_i) be the length n prefix of the sequence 
r. For any given multi-sequences r = (ri,-- - ,rj,--- ,rm), where each 
''i ~ {cj.t}t>o is a sequence over F, we identify it with the element r = 
(ri(z), r2(z), • • • ,rm{z)Y G F{{z~^))^ with valuation larger than 0, where 
rj{z) = ^t>o^j,t'^~^~^ i^ formal Laurent series identified with the jth 
sequence rj, and let A"'^ = {r^\r^\--- ,rm^Y be the length n prefix 
of the multi-sequences r. Given a polynomial q{z) over F, we call it a 
characteristic polynomial of A"'^ if Iv{r — -^^^) > {nT-,n), or equivalently 

Iv{{rq{z)}) > {m,n — deg{q{z))); and call it a minimal polynomial of r^") 
if it is a characteristic polynomial of of the smallest degree; and call 
deg{q{z)) the linear complexity of r^"), denoted by Ln{r), if q{z) is a min- 
imal polynomial of r^"). The multi-sequence synthesis problem is: Given a 
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multi-sequences r € F{{z~^))"^ , find a minimal polynomial and the linear 
complexity of r^"^ for each n > 1. 

This problem is solved by using m-continued fractions as shown in the 
following corollary. 

Corollary 6.11. Given a multi-sequences r E F{{z~^))"^ , let C E T(r). 
Write C in the form of (jS-lj) . Then qk{z) is a minimal polynomial of r^") 
and Ln{r) = dk for Uk < n < Uk+i, < k < to, where we let = oo in the 
case uj < oo, and no = 1. 

Proof. It is easy to prove. □ 
7. Proofs of Main Results 
In this section we complete the proofs for all results of this paper. 
7.1. Approximants of m-continued fraction-II. . 

Proof of Theorem 16.11 : Keep all the notations made for C. We prove it 
by induction on k. For k = 1, we have qi{z) = ai^hj^{z) by Proposition 15.31 
hence deg{qi{z)) = ti = di. Assume we are done for < k, where k > 2. 
From Proposition 15.31 we have 

qkiz) = qk-i{z)ak,hk{z) + Qk-2,hk + Qk-2,jak,j{z). 

j¥'hk,l<j<m 

The wanted result deg{qk{z)) = d^ follows by observing the following facts: 

• deg{qk^i{z)ak,hk{^)) = dk-i + tk = dk (induction assumption). 

• If Qk-2,j 7^ 0, from Proposition 15.31 we see l{k — > 1, then 
deg{Qk-2,j) = deg(g,(fe_ij)_i) = ( by induction assump- 
tion). In particular, deg{Qk-2M) < dk- 

• For j 7^ hk and Qk-2,j 7^ 0) from the above and Proposition 15.61 we 
getdeg{Qk-2,jak,j{z)) < then deg(Qfc_2jafcj(z)) < 
dk, since l{k — l,j) = l{k,j) in this case. □ 

Lemma 7.1. Let C, which is defined as J5. ij) . be an m-continued fraction, 
and r = 93(C), and keep all notations made for it. Let 

V^{z) P^{z) 
qr{z) qt{z) 
lr,t = P^{z)qt{z) - p^{z)qr{z), yO<r<t. 

Then 

(1) /w(z,. i) = {h,v), if and only if Lv{y^^^) = {h,v + dr + dt). 

(2) The following statements are equivalent to each other: 

^^(Ht-iJ = (^t' ifO<t<k, 

I^iz^^i^t) = iht,vt - dt), if <t < k, 

dv{y^.^) = {hr+i,nr+i), if <r <t < k, 

Ivizr^t) = ihr+i,Vr+i - dt), if < r < t < k. 
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Proof. The item (1) can be verified easily by Theorem IH. 11 The item (2) can 
be obtained by a routing proof. □ 

Proof of Theorem 16.21 : Keep aU the notations made for C. It is enough 
to prove that Iv{y^_^^ = {ht,nt) for t>l, which is equivalent to 

(7.1) Iv{zt_,^t) = {ht,vt-dt), l<t 

from Lemma [7.11 For t = 1, Zq i = —p^{z), then Iv{zq i) = Iv{p^{z)) = 
{hi,di — ti) = {hi,vi — di) by Corollarv lb.31 Suppose we are done for t < k, 
where k >2, which together with Lemma l7 . 1 1 implies 

(7.2) Iv{zr,t) = {hr+l,Vr+l - dt), < r < t < k. 
By Proposition 15.31 we have 

= Vf^_^{'ik^iak,hk{^) + Qk-2,jak,j{z) +Qk-2,hk) 

-iPk-l"'k,hki^) + Pk-2,j"'k,j{z) + Pk^2Mlk-l 

= -{Pk-2,hk^k~l - Pk-iQk-2,hk) 

- {Pk-2,j(lk-l - P^_^Qk-2,j)ak,j{z) 

= ~^k,hk ~ ^ Z_i^jakj{z), 

where 

Z.k,j = Pk-2,j(lk-i -Pf^^^Qk-2,j VI < j < m. 
By Proposition 15.31 we can get 

if lik-lj)>l, 
ejqk-i{z) if l{k-l,j)=0. 

Hence by (|7.2|1 we get 

(7.4) Iv{Zkj) = (j, vi(^k-i,j) - dk-i). 
In particular, we have 

(7.5) Iv{Z_k,hJ = {hk,vi(k~i,hk) - dk-i) = {hk,Vk - dk). 
For j 7^ /ifc, we have 

Iv{Z^jak,j{z)) > {j,vnk-i,j) - dk-i - {v^kj) - vi(k-i,hk) ~ ^hk,j)) 

= (i) Vi(^k-l,hk) ~ ^k-i + Ihkj) 

> {hk-,Vl(k-l,hk) ~ (^k-i) 

= ihk,Vk-dk), 



^k,j 
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where the first hne comes from (|7.4j) and Proposition l5.(il (vi(!-j) = ^fc,j)> the 
second line comes from l{k,j) = l{k — l,j) since j 7^ h^, the last two lines 
come from definitions. This together with (|7.3j) and (|7.5j) leads to equation 
((ZU when t = k. □ 

Lemma 7.2. Let C, which is defined as i5. 1]) . be an m- continued fraction, 
and r = (p{C), and keep all notations made for it. Assume Ci{z) E F[z\, 
deg(ci(2;)) < tj+i, < « < u; — 1. Then 

(1) Iv{{rqi{z)ci{z)}) = {hi+i,Vi+i - deg(cj(z))). 

(2) Iv{{rqi{z)ci{z)}) 7^ Iv{{rqj{z)cj{z)}) for all < j ^ i < uj and 
Ci{z)cj{z) 0. 

Proof. It is easy to prove. □ 

Proof of Theorem l6.4t Keep all the notations made for C. This theorem 
is true in the case k + 1 = uj < 00, since r = \ in this case. We need only 
consider the case A; + 1 < a;. Denote d = Aeg{q[z)). We have dk < d < dk+i 
by the assumption. It is enough to prove Iv{{rq{z)}) < {hk+i,nk+i — d). 
We can write q{z) = Z]o<i<fc where (leg{ci{z)) < deg(gi+i(z)) - 

deg{qi{z)) = tj+i for i > (note that qo{z) = 1) and deg(cfc(z)) = d—d^ > 0. 

Note that {rq{z)} = Yl {LQiiz)ci{z)}, from Lemma |7.2I and Lemma \'A.^\ 

0<i<k 

we have 

Iv{{rq{z)}) = min{Iv{{rqi{z)ci{z)}) | Ci{z) ^0,0 <i<k} 
< Iv{{rqk{z)ck{z)}) 

= {hk+i,Vk+i - deg(cfc(z))) = {hk+i,nk+i - d).n 

7.2. m-CF Transform and m-continued fractions-II. In this subsec- 
tion, at first we restate the m-CF Transform in terms of matrices, then we 
prove Theorem I6.1UI Before giving the matrix version of m-CF Transform 
we give some necessary definitions. 

Definition 7.3. {Base matrix and its D-component ) 

Let R he a matrix of order m over F{{z~^)), we call it a base matrix if 
each of the columns of R is nonzero, and the index of the j-th column of 
R is J for each 1 < i < m. For a base matrix R, the D-matrix A = 
Diag.{z~'^^ , z~^^, ■ ■ ■ , z^'^™) will be called the D-component of R, where cj 
denotes the valuation of the j-th column of R. 

Proposition 7.4. (Matrix version of m-CF Transform) The m-CF 
Transform which is introduced in subsection 6.2. can be restated in terms 
of matrices as follows: 

Given r G F{{z~^))"^, r / and v{r) > 0. Initially let R-i = Im, Lo = L- 
Suppose {—Rk-2 Lk-i) ^ -^m,m+i(-^((-z > 1, is defined, where 

Rk-2 G M„,^(F((z-i))) andr^Ji G F{{z-^)Y\ Iflfc.i =0, the algorithm 
terminates and denote uj = k. If r^_i 7^ 0, then do the following steps: 
(1) Let hk = I{rk-i) G Z^- 
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(2) Let Rk-1 be the matrix of order m and Sfc„2 ^ F{{z~^))'^ such that 
{-Rk-i Sfc_2) = {-Rk-2 r,^-i)Ehk G Mm,ni+i{F{{z-^))), where 

is defined in Notation 15.21 

(3) Let i?fc^iSfc_2 = Ak + ak G F((z-^))™, where ^fc = L^fc^iSfc-2J' and 

afc = {Rk-iSk-2}- 

(4) Let be any one in such that a^, = Ak — e^, where ej^ € 
/^(Afc^iafc) < /^(Afe-iefc), and Afc_i is the D-component of Rk-i- 

(5) Let G such that = ife_2)^(afe) e 

Denote w = oo if the above procedure never terminates. 

Proof. It follows directly from the following two lemmas (Lemma 17.61 and 
Lemma 17.51 ) . □ 

Lemma 7.5. Keep notations made in Proposition \ 7.4\ Denote Pl_i = 
ak-i + efc-i. For k > 1, we have 

(1) Pk-i = ^/ ™^ '^''^^y if Lk~i — 0- 

(2) Rk-i is a base matrix. 

(3) r,_i = i?fe_2/3Li- 

(4) v{^k-i,) = [''\^'~'i , % 

^ ^ ^ '^^ \ v(Afc_2/3fc_i) J = hk- 

(5) Denote R^\s,_, = (• • • ,p^,,,---r G Fi{z-')r, p^. G 

T/ien pl - = i/ j / hk, pIj^^ = -pi^^, where (3l_^^- is the 

j-th component of (5^_^. 

Proof. It is easy to check. □ 

Lemma 7.6. Any base matrix R is invertible, and Iv{Rr) = Iv{Ar) for all 
r G F[{z~^))^ , where A is the D-component of R. 

Proof. It is easy to prove. □ 

We divide Theorem I6.1UI in the following two parts, then prove them 
separately. 

Part 1 of Theorem lHTTUl T(r) C C(r), and i/;(r) is a strict m-continued 
fraction. 

Part 2 of Theorem IH7TU1 C(r) C T(r), Moreover, if C G C(r) is strict, 
then C = ip{r). 

To start the proof, we make a notation which will be used later. 
Notation 7.7. Let i? be a base matrix, define 
(7.6) v{R) = (ci,C2,--- ,Cj,--- ,Cm), 

where cj denotes the valuation of the j-th column of R. We use the notation 
{h,v) > v{R) to denote the fact that v > c^. Notations {h,v) < v{R), 
v{R) > {h,v), etc. should be interpreted similarly. 
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Lemma 7.8. Let A be a D-matrix, / r G Then Iv{Ar) > 

v{A) ifv{r) > ( including r = 0); and Iv{Ar) < v{A) if / r G F[z]"^. 

Proof. It is easy to prove. □ 

Lemma 7.9. Let R be a base matrix, A be the D-component of R, s £ 
Fiiz-^))"", A = [R-^s\, and a = {R-^s}. Let e (£ Ffz]"^, anda = A-e. 
Then the following conditions are equivalent: 

(1) Iv{Aa) < Iv{Ae) . 

(2) a = e = 0, or Lv{Aa) < Iv{Ae). 

(3) Iv{-Ra + s) = Lv{Aa). 

(4) Iv{-Ra + s) > v{A) 

Proof. It is easy to prove. □ 

Lemma 7.10. Let R be a base matrix, A be the D-component of R, s G 
F{{z-^))'^, A= [R-^s\, and a = {R-h}. Lete(^F[zf^, and a = A-e. 
Ifv{R) > Iv{s) and Lv{Aa) < Lv{Ae), then 

(1) iv{s) = Lv{Aa) = Lv{AA), and a / 0. 

(2) Iv{Aa) < Lv{-Ra + s) > v{A). 

(3) a = if and only if Supp^(Aa) < Iv{—Ra + s) > v{A). 

Proof. It is easy to prove. □ 

Lemma 7.11. Let r G F{{z~^))^ , v{r) > 0, and keep the quantities ob- 
tained by acting the m-CF transform on r, say, hk, Rk-i, ^fc-2> '^k, '^k-i, 
^k, QLk; Lk ^'^'^ 6tc. Let 

C = [ao,hi,a^,h2,a2, - ■ ■ ,hk,ak,---] G T(r). 
Then, for k : 1 < k < to we have 

(1) Rk-i is a base matrix, so the m-CF transform is well defined. 

(2) v{Rk-i) > Lv{sk^2) = Lv{Ak-iak) < Lv{r,,) > v{Ak-i). 

(3) deg{ak,hk{^)) ^ 1; ^-fi-, C satisfies the condition 1. 

(4) Based on the item (3), we may keep the quantities associated with 
C , which are defined as in section 5, say, tk, Dk, Vkj, Vk, etc. Then 

Iv{Ak-ij) = Iv{Rk-i,j) = ij,Vk,j)- 
Afc_i = Dk- 
Iv{u^i) = {hk,Vk). 

Iv{Sk-2) = ihk,Vk-l,hk). 

Proof. It is easy to prove. □ 
Proof of Part 1 of Theorem IHIIJI 

From Lemma fT.lll whose notations are kept, we see C satisfies the conditions 
1, so, in order to prove C is an m-continued fraction, it is enough to prove 
C satisfies the conditions 2 and 3. Prom Lemma 17.111 we can see: 

• For A; : 1 < /c < u; - 1, we have (/ifc, ■Wfc-i.hJ = Lv{s^_2) < Lv{r^) = 
(/ifc+i, ffc+i) , which means condition 2 is satisfied. 
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• For /c : 1 < A; < cj, we have (/ifc, Wfc-i,/ife) = Iv{sf,_2) = If (Afc-ia^) = 
Iv{Dkaj^) , which shows that condition 3 is satisfied. 

In order to prove C is strict if C = ip{r), it is enough to prove C satisfies 
the condition 4. From C = ^(r) we have = \_Rk-i-k-2\A,, i- From 
Lemma l7. Ill we see Rk-i is a base matrix, and v{Rk-i) > Iv{s_^_2)', from 
definition of m-CF transform we see Iv{llk-i{R^\sj^_2\) < Iv{Ak_iej^), 
so we can apply Lemma [7. 101 ( taking R = Rk-i,s = Sfe_2)i = ffc ) to get 
Supp^iAk^iOk) < /w(--Rfc-iafc + Sfc_2) = Ivizf^). Then, by LemmaEIU 
we get S'npp+(Afc_iafc) < Iv{r_f^) = (/ifc+i, ti^+i) for 1 < A; < lj — 1, i.e., the 
condition 4 is satisfied. 

Finahy we prove V3(C) = r. Keep all the notations made for C in 
section 5. By definition of m-CF transform, we can see {—Rk-i,Zk) — 
{-Rk-2,Lk-i)Eyik^{ak) = ■■■ = {-Im,r)Bk. This leads to = -p^ +mk 
(the last column of Bk is {Pj^,qkY), so r— p^/qk = Lk/lk- By Lemma IV. Ill 
we have Iv{ri,) = {hk+i,Vk+i), then Iv{r-pjqk) = (/ifc+i, ^'fc+l + deg(gfe) = 
{hk+i,Vk+i + dk). Therefore Ximk^ooPf^/ qk = L since Imik^oodk = 00, in 
other words, ip{C) = r. □ 

Lemma 7.12. Let C = [qq) ^Ij ^i) ^2, «2' ' ' ' ■> ^k-,(kk-, ' ■ ■ ]; which is defined as 
i5.1]) . be an m-continued fraction, and keep all notations made for it, say, Vk, 
Dk, etc. Letr = ip{C). Define Rk-i G M^,^(F((z-i))), r^, G 
and Sk^2 ^ F({z~'^))"^ iteratively as below ( 1 < k < uj ): 

{R-i = Im, Lo=L 
{-Rk-1 Sk,2) = {-Rk-2 rj,_,)Eh, e M^^m+i{F{{z-^))) 
{-Rk-1 rj,) = {-Rk-i Sk_2)A{ak)&Mm,m+i{F{{z-^))). 

Denote by Rk-ij thej-th column of Rk~i, denote by Ak-i the D-component 
of Rk~i, and denote Ck = lRk\sk-2\ Then, for 1 < k < uj, we have 

(1) IviLk-i) = {hk,Vk). 

(2) Rk^i is a base matrix, and Ak-i = Dk- 

(3) Ivizk) = Iv{-Rk-iaj, + Sk^2) > v{^k-i)- 

(4) v{Rk-i) > Iv{sk_2)- 

(5) Iv{Ak-i{R^\sk_2})<Iv{Ak-ie^). 

Proof. It is easy to prove. □ 
Proof of Part 2 of Theorem IHUJl 

Keep notations made in Lemma 17.121 Prom Lemma 17.121 we have /(r^_i) = 
hk and Iv{Ak^i{R]^\s_k_2}) < Ivi^k-i§.k) for 1 < < w, which show that 
1)7. 7() is an m-CF transform, so C € T(r). To prove C = ip{r) when C is 
strict, it is enough to show Of, = L-^fe-ilfc-2J Aj. ^ for A; : 1 < A; < lj. In fact, 
if C is strict, we have Supp'^{Dkak) < (^fc+i) ^fc+i) for 1 < A; < u; — 1, then 



Supp'^ (Ak-ittk) < {hk+i,Vk+i) = Ivir^) = Iv{-Rk-iak + Sk_2) > v{Ak-i) 



MULTIDIMENSIONAL CONTINUED FRACTION AND RATIONAL APPROXIMATIONL7 



by Lemma l7.12| thus Uj^ = lRk-i^k-2\Ak i according to Lemma l7.1()l ( 
taking R = Rk-i,s = Sj.-2-,§. = ffc )) since all the conditions for applying 
Lemma I7.1UI are satisfied by Lemma 17.121 □ 
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